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Abstract 

New FOCUS data on D + — > K~ir + ir + are fitted to the k, together with earlier data 
from LASS, E791 and BES 2. There is a clear low mass Kir peak due to the k pole. 
Inclusion of the 7 = 3/2 Kir amplitude gives only a marginal improvement to the fit and 
almost no change to the K peak. An improved formula for the K gives a better fit than 
that used earlier. The k pole moves to 663 ± 8(stat) ± 34(syst) — i[329 ± 5(stat) ± 22(syst)] 
MeV. The Kq(1430) pole is at 1427 ± A(stat) ± I3(stat) - i [135 ± 5{stat) ± 20(syst)} MeV. 

PACS: 14.40.-n, 13.75.Lb, 11.30.Qc 
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1 Introduction 

There are two objectives in this paper. The first is to present a fit to new FOCUS data [1] and 
compare them with earlier E791 data [2]. There are some systematic discrepancies between these 
two sets of data, although their effects are minor. Like E791, the FOCUS group has determined 
the K~ir + S-wave amplitude in magnitude and phase in 40 mass bins covering the whole mass 
range available in D + — > K~7t + tt + . These data are fitted here. 

The second objective is to introduce a necessary modification to the formula used earlier to fit 
the k. Both a and k resonances have widths which are strongly s-dependent, giving them unusual 
features. This s-dependence originates directly from chiral symmetry breaking. Hopefully, the 
formulae developed here will be useful to experimental groups as direct replacements for the 
usual Breit-Wigner amplitudes used to fit J p = + resonances such as k and A^ (1430) in 
Dalitz plots. The parametrisation for the Kir channel is highly convergent and well determined 
in terms of just three parameters, one of them related closely to f%. Some allowance is also 
required for coupling to Kr] and Kr]', but in the absence of direct experimental information 
on these channels, this requires just one coupling constant for each channel. If data for these 
channels become available, the present treatment of Kir could be substituted also for Kr] and 
Kr] 1 . 



2 Formulae 

The key point about both a and k is the Adler zero in the elastic amplitude at s ~ rn 2 K — for 
the k and at s ~ 7^4/2 for the a. These zeros are crucial features of chiral symmetry breaking. 
The mechanism of this symmetry breaking is fully understood [3] [1] [5] and is confirmed in 
outline by Lattice QCD calculations [Bj. References [3-5] give illuminating discussions of the 
detailed mechanism. Pioneering work on the k and a was done by Pelaez, Oset and Oiler and 
collaborators [7] [S] [S] HDJ. 
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The consequence of the Adler zero is that the Ktt amplitude rises nearly linearly with s near 
threshold. In Ref. [IT] , the k amplitude was parametrised as a Breit-Wigner amplitude with 
s-dependent width: 

N(s) MT el 
Jel D(s) M 2 -s- iMT total 1 ' 

MT ei = exp(-ak 2 )g 2 (Kn)p Kn (s); (2) 

M z — s A 

here s A = 0.2367 GeV 2 is the mean position of the Adler zero for K + tt~ and K°7T°, and p(s) is 
Lorentz invariant phase space 2k/ y^, where k is centre of mass momentum. Also g are coupling 
constants and a is a fitted constant. The exponential form factor in (2) accommodates the 
experimental fact that T e i gradually flattens off at large s. 

There are two weaknesses in these formulae. The important one is that they assume the k 
phase shift eventually reaches 90°, although it was found in [11] that this was only at ~ 3.3 
GeV, well above the mass range of available data. The FOCUS data, combined with earlier 
data, now reveal that the k phase shift appears to reach only ~ 55° at 1.5 GeV and may never 
reach 90°. It is straightforward to modify the formulae to cater for this possibility. The second 
weakness is that the mass and width of the pole itself are far removed from M and T to tai of Eqs. 
(1) and (2), for reasons described in Section 4. Formulae can be rewritten to be more closely 
related to the k pole itself. 

A better form of the equations may be obtained by dividing both numerator and denominator 
of eqn. (1) by M 2 and writing 

, = Mf ~ s A )F 1 (s)p K7[ (s) 

hl 1 - s/M* - i Ej h{s - 8^(3)^(8) ■ 1 } 

For the k, the term s/M 2 is small; bj are constants and the summation in the denominator is 
over j — 1 — 3 for Kn, Krj and Krf channels in principle, though the coupling to Kt] turns out 
to be insignificant. Form factors Fj are discussed in detail below. 

Weinberg predicted the scattering lengths of pions from any target [12]; this sets the scale of 
chiral symmetry breaking, hence b\ in terms of f 2 . Eq. (3) can be recast so as to expose the 
scattering length explicitly. This requires a transformation in s of the denominator of (3). The 
algebraic manipulations are to write 1 — s/M 2 = 1 — As ~ [1 — A(s — s^ r )]/(l + Asthr) for small 
A, next multiply top and bottom of (3) by (1 + As t hr) then replace all bj by Bj = bj(l + As t hr)] 
here s t hr are evaluated at the threshold of each channel j. The result is 

_ N(s) Bx(s - s A )F d p Kn (s) 

Ul ~ D(s) ~ 1 - A(s - s thr ) - % Ej B,{s - s Aj )Fj(s) Pj (s) ' 1 ' 

If A is positive, Re D eventually goes to zero and the elastic phase shift reaches 90°. The case 
where A is negative accomodates the possibility that the phase shift does not reach 90°. The 
denominator however still allows a k pole. The scattering length a is 

a = 2B 1 (s t hr ~ S A )/Vsthr- ( 5 ) 

There is an alternative way of viewing this formula. The factor 1 — A(s — sthr) in the 
denominator of Eq. (4) may be regarded as an empirical form which exhibits the scattering 
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length explicitly and parametrises successfully the real part of the amplitude in a region close to 
threshold and the Adler zero. For elastic scattering, a form factor allowing a controlled departure 
of the numerator N(s) from linearity is F± = exp(-aikf), where a% is a fitted parameter. Note 
that the exponential dependence of the form factor combines both a conventional form factor 
for the radius of interaction and an empirical departure of the numerator from linearity; values 
of ot\ and A are accurately determined by FOCUS and E791 data. 

The Kr] channel turns out to be negligible; for Kr]', the Adler zero is at s = m 2 , —m 2 K /2. The 
form factor for Krj' is not known, because of lack of data for this channel. Above its threshold, 
the value <x,- = 4.5 (GeV/c)~ 2 is adopted from a wide range of Crystal Barrel and BES II data; 
this value corresponds to a reasonable radius of interaction 0.72 fm. Above the Krj' threshold, 
small changes in ctj may be taken up by small alterations to parameters fitted to -fTo(1430) and 
-fTo(1950) and the fit to data changes by less than the errors. 

Below the inelastic threshold, the Flatte prescription is adopted, continuing p analytically: 
P — > i\p\ [13] ■ However, \p\ increases below threshold and requires a cut-off; otherwise, a myriad 
of open channels at high mass dominate Kir elastic scattering at low mass. Such a prescription 
would obviously be unphysical. With the accuracy of present data, any reasonable cut-off will 
do, and the simple one 

Fj = exp[-oy|A;?|] (6) 

fits data adequately. The best source of information on sub-threshold form factors comes from 
Kloe data on — > 77r°7r°, where a similar cut-off is required for the sub-threshold a — > KK 
amplitude [TJ]. Rather a sharp cut-off is required for those data and optimises with atj = 8.4 
(GeV/c) -2 , though a value as low as half this is acceptable. For present data, a similar form 
factor is definitely required. If a value ay < 2 (GeV/c)~ 2 is used, the ifo(1430) — > Krj' amplitude 
near the Ktx elastic threshold becomes unreasonably large. There is a weak optimum at oij = 4.5 
(GeV/c) -2 . With this value, the sub-threshold contribution rises rapidly from the Krj' threshold 
to a Kn mass of 1400 MeV and thereafter varies slowly and smoothly down to the K-k threshold. 
Any alteration to this value of ctj is absorbed by small changes to the fitted parameters of 
i£o(1430) and small changes to fitted values of a± and Bj. The extrapolation of the amplitude 
to the pole is stable: the real and imaginary parts of the pole position change by 20-30 MeV for 
variations of <Xj in the range 3.5-6.5 (GeV/c) -2 . This will be taken as a systematic error on the 
pole position. [If data one day become available on the form factor for Krj' above threshold, it 
will be possible to calculate the sub-threshold continuation from a dispersion relation like Eq. 
(8) discussed below, but presently this is academic] 

For production processes such as D — > (Ktt)ii, the denominator D(s) of the amplitude must 
be identical to that of elastic scattering by Watson's theorem [15J. However, the numerator 
can be (and is) very different. Data on J/^> — > uj-k + t[~ [16] and K + -k~K~it + [17] [18] and 
present data on D — > (Kit)?! may be fitted accurately taking N(s) as a constant. The result 
is that cr and k poles appear clearly in production data, but are obscured in elastic scattering 
by the nearby Adler zero in N e i(s). Precise theoretical work on elastic scattering using the Roy 
equations does however reveal a and k poles clearly [19] [20J. 

The Adler zero is a property of the full Kw elastic amplitude and therefore needs to be 
included into the widths of iT (1430) and K o (1950). This is done using Eqs. (1) and (2) for 
each of iT (1430) and K Q (1950); the exponential form factor is taken to be the same as for the k. 
Although Z£o(1430) may be fitted with a Breit-Wigner amplitude of constant width, the resulting 



3 



contribution to the scattering length is far too large to agree with Weinberg's prediction. 

In order to satisfy unitarity for LASS data, the full S-matrices of k, Xo(1430), i£o(1950) and 
the Kir 1 = 3/2 amplitude are multiplied. This prescription is not unique, but gives a good fit 
to the interference between k and K (1A30). 

The FOCUS, E791 and BES 2 groups have all determined the magnitude and phase of the 
Ktc S-wave with respect to other strong components in the data. To fit production data, the 
isobar model is used, with a complex coupling constant in the numerator of each amplitude 
instead of MT e i, specifically Aexp(z0)/(M 2 — s — iMT tot ). The first two sets of data require a 
phase difference between k and A" (1430) ~ 75° different to elastic scattering. 

One alternative procedure for fitting the data has been explored, but in practice turns out 
to give no better result than the Flatte prescription, though it does provide a cross-check. In 
principle, more complete formulae for a Breit-Wigner resonance are [21] 



f = MTa 

Jd M*-s-m{s)-iMT total {sy K) 
M 2 — s f ds' T total (s') 

The scattering amplitude is an analytic function of s. Any s-dependence in MT to tai gives rise to 
a corresponding contribution to the real part of the denominator via m(s). The rapid opening of 
the Kr]' threshold generates a sharp spike in m(s) precisely at the Kr]' threshold. A subtraction 
in (8) at mass M improves the convergence of the integral. 

The problem which arises is that this spike depends on the mass resolution of the experiments. 
Without an accurate knowledge of mass resolution and possible variations over the Dalitz plot, 
inclusion of m(s) is impractical for Kr] and Kr]' channels, though satisfactory for the Ktt channel. 
The Flatte formula only partially reproduces the effect of m(s). However, in practice, it allows 
small adjustments of M and gjc^/gjc^ which fudge the effect of m(s) by optimising the fit to 
data. Readers interested in the effect of m(s) are referred to Ref. [22J, where it has been included 
into fits to ao(980) in Crystal Barrel data. In that case, the mass resolution is known accurately 
and included. It smears out the spike at the KK threshold seriously, as shown in Fig. 6(b) of 
that paper. Even then, the dispersion integral for m(s) is sensitive to the form factor, and the 
height of the spike needs to be fitted empirically. Here, the dispersive correction will be used 
for the Kir channel so as to assess possible systematic errors in parametrising that channel. 



3 The fit to data 

An important feature of present work is that data from LASS for elastic scattering [23], FOCUS 
[T] and E791 [2] and BES 2 [TS] are fitted simultaneously with consistent parameters. This 
reveals the strengths and weaknesses of each individual set of data. In BES 2 data, there is a 
strong i^o(1430) peak, together with a known smaller contribution from the narrower 7^2(1430). 
This peak is included in the present fit since it gives the best determination of parameters for 
AV1430). 

The E791 data were published including a form factor in N(s) for production, but it was 
subsequently shown [21] that the fit optimises with this form factor set to 1. That correction 
has been applied to data fitted here. FOCUS also take the form factor for production to be 
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Figure 1: The result of the combined fit for (a) magnitudes and (b) phases of FOCUS amplitudes; 
(c) and (d) the corresponding fit to E791 data 



1. Figs. l(a)-(d) show the fit to FOCUS and E791 data. One sees some modest systematic 
discrepancies between them. The fit to either may be improved by omitting the other, but both 
are included in the final fit. Values of x 2 will be given for a variety of combinations in Table 
2; x 2 is calculated directly from FOCUS and E791 tabulations, combining their statistical and 
systematic errors in quadrature. 

Figs. 2(a) and (b) show the fit to LASS magnitudes and phases. It is well known that these 
data stray slightly outside the unitarity circle above 1.5 GeV, so errors in this region have been 
increased until an average x 2 of 1 is achieved; this refinement has no significant effect on the fit. 
Figs. 2(c) and (d) show the fit to BES 2 data. The phases of two high mass points now look 
high; however, since they are direct measurements of phase in the data, they are retained in the 
fit. The fit to the 1430 MeV peak in BES 2 data is shown in Fig. 3(a). 

The fit to all data requires small contributions from K Q (1950). Its parameters have been 
allowed to vary by la from averages quoted by the Particle Data Group [23|. They finish at 
M = 1967.4 MeV, Tel = 115 MeV, T tota i = 376 MeV, but these changes have insignificant effect 
on parameters fitted to k and i£o(1450). 

There is one quite large change in the fit to ifo(1430) compared to that reported in Ref. 
[24] . The ratio gxri'/dKn reported there was 1.15. This now falls to 0.62 ± 0.06. This change 
may be traced directly to the improvement in Eqs. (3-5) over Eqs. (1-2), used earlier. The 
latter accommodated a Kir phase of ~ 72° in BES 2 data, but forced a large value of g\ , for 
_Ko(1430) via the sub-threshold contribution of the Krj' channel. However, the earlier fit was 
visibly not perfect from 1.15 to 1.35 GeV because this contribution was too large. That problem 



5 




0.75 1 1.25 1.5 1.75 



c 1 

S) 
ro 
E 

0.75 
0.5 
0.25 





3ES 2 



_J i I i I i I i L 

0.75 1 1.25 1.5 1.75 
Mass Kit (GeV) 



150 



50 - 




0.75 1 1.25 1.5 1.75 



40 











IT l I 


, I , I , I 



0.75 1 1.25 1.5 1.75 
Mass Kit (GeV) 



Figure 2: The result of the combined fit for (a) magnitudes and (b) phases of LASS amplitudes; 
(c) and (d) fits to BES 2 data. 



now disappears. With the new formulae, the phase shift of the k peaks at 55° at 1530 MeV, 
though any fall between there and the end of the mass range is within statistics. The magnitude 
and phase of the k contribution are displayed in Figs. 3(b) and (c). The vertical scale for Fig. 
3(b) is arbitrary. 

Inclusion of the coupling of k to Kr] improves x 2 by only 4.6, i.e. 2a. It also tends to 
destabilise the fit because of interferences between the sub-threshold Kr] component and Kir. It 
is therefore omitted from the final fit. There is now evidence for a component due to k — > Kr]' . 
It optimises at g 2 {r]' K) / g 2 (K-n) = 0.085 ± 0.020 and improves the fit by 14.6, i.e 3.8er. It is not 
strongly correlated with the contribution from k — > Kr]' . The small structure in the k phase 
near 1.5 GeV is due to coupling to Kr]'. Direct data on Kr]' would be valuable to improve the 
determination of the coupling of k and iT (1430) to this channel. 

3.1 Possible effects of the Kit 1 = 3/2 amplitude 

This amplitude has been parametrised by Pennington. His formula for the K-matrix, quoted in 
Ref. |2S], will be repeated here for completeness: 

K 3/2 = ( g ~°- 27 ) (-0,22147 + 0.026637s - 0.00092057s 2 ), (9) 

° norm 

where s = s/s norm — 1, s norm = {m 2 K + m 2 ,) and s is in GeV 2 . Inclusion of this amplitude 
improves x 2 only by 7.7. This is an indication that the data agree well with both the magnitude 
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Figure 3: The fit to (a) the intensity of the 1430 MeV peak in BES 2 data; (b) intensity and (c) 
phase of the fitted n amplitude. In (c), the dashed curve shows the shape of the LASS effective 
range formula. 



and the phase of the broad k amplitude. The improvement is not concentrated in any particular 
mass region or data set, but is fragmented into small improvements distributed almost randomly. 
What happens in the fit is that destructive interference develops between k and the 1 = 3/2 
amplitude. These are familar symptoms that two broad amplitudes are combining to fit minor 
defects in the data. 

There is a physics reason why the contribution from repulsive amplitudes should be small. For 
attractive interactions (particularly resonances), the wave function is sucked into small radii, r; 
for repulsive interactions, it is repelled to large r by the potential barrier and therefore reduced 
in magnitude. There is evidence that J/ty and D decays involve short-range interactions. This 
comes from the absence of form factors in their decay processes. From the earlier analysis of 
E791 data, the interaction was found to have an RMS radius < 0.38 fm with 95% confidence. 

In an earlier FOCUS publication on D + — > K~tt + 7t + , the repulsive 1 = 3/2 Kit amplitude was 
included, but led to massive destructive interference with the k. In the earlier FOCUS analysis, 
the Adler zero was not included into the k amplitude, which was parametrised by a broad 
Breit-Wigner resonance and a constant interfering background. With three broad amplitudes, 
destructive interferences are able to patch up minor defects all over the Dalitz plot. Data on 
other charge combinations could identify production of I = 3/2 Kir and I = 2 tttt amplitudes. 

The Cleo C collaboration has also presented an analysis of D + — > K~7T + 7r + where they include 
the I = 2 7T7T amplitude [27]. They present three solutions, all of which contain huge 1 = 2 
amplitudes with an intensity a factor ~ 30 larger than K*(980)tt. This is clearly symptomatic of 
large destructive interferences with the k amplitude, since the if* (890) is clearly visible by eye 
in the Dalitz plot and Kir mass projection. Their second solution is made with an amplitude 
derived directly from a complex pole at M = 706.0±1.8(stat)±22.8(syst)-i(319A±2.2(stat)± 
20.2(syst) MeV. In view of the huge interference with the 1 = 2 amplitude, it is difficult to know 
how reliable this result is. 
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Parameter 


Value 


K 


A 
si 


— U.UoU it u.uoz 




D 
D\ 


Z.O/o it U.Uoy 




OL\ 


U.000 ± U.Uoo 




2 12 
9 Kn'l 9 K-k 




^o(1430) 


M 


1.479 ±0.004 




9Kit 


0.284 ±0.012 




9kt)'I 9 Kit 


0.62 ±0.06 



Table 1: Parameters fitted to k and iTo(1430) in units of GeV 
3.2 Pole positions 

Table 1 gives fitted parameters for k and _Ko(1430) in units of GeV. Pole positions are remarkably 
stable if individual sets of data are dropped from the fit, or even pairs of sets. As an illustration, 
Table 2 shows values when each of the sets of data listed in column 1 are dropped. There are 
14 fitted parameters: the 7 listed in Table 1, 6 for three complex coupling constants fitted to k, 
_Ko(1430) and fTo(1950) in D — * Knir, and one scale factor for the magnitude fitted to i^o(1430) 
in BES 2 data. The total x 2 is 349.4 for 272 degrees of freedom. Errors quoted in Table 1 
for each parameter correspond to changes in x 2 of 1.28, the mean value per degree of freedom. 
However, there are some correlations between parameters, with the result that the quoted errors 
overestimate somewhat the changes expected in Table 2. The final two columns in Table 2 show 
changes in \ 2 when each set of data is dropped and also the number of points dropped. 



Dropped 




K 


K (1430) 


Ax 2 


points 


LASS phases 


661 


-i325 


1432-H27 


34.9 


41 


LASS magnitudes 


663 


-1334 


1435-H40 


44.9 


41 


FOCUS magnitudes 


663 


-i329 


1427-H37 


41.7 


40 


FOCUS phases 


660 


-i327 


1427-H29 


98.2 


40 


E791 magnitudes 


661 


-i328 


1426-H35 


73.2 


38 


E791 phases 


663 


-i329 


1424-il39 


48.2 


38 


BES K (1430) 


675 


-i341 


1431-H44 


54.0 


29 


BES magnitudes 


640 


-i321 


1425-il32 


21.3 


9 


BES phases 


657 


-i323 


1426-il33 


27.4 


10 


FOCUS all 


660 


-i328 


1425-H32 


158.4 


80 


E791 all 


662 


-1327 


1424-H40 


126.4 


76 



Table 2: Pole positions in MeV when individual sets of data are dropped; the last two columns 
show the change in x 2 of the fit and the number of points dropped. 

Using all sets of data, the i^ (1430) pole position is 

M = 1427 ± A(stat) ± 13(syst) - z[135 ± 5(stat) ± 20(syst)] MeV. (10) 
The systematic error is compounded from the worst case and from uncertainties in relative 
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contributions of ^(1430) and A" (1430) (and their relative phase) to BES 2 data on the 1430 
MeV peak. For the k, the pole position is 

M = 663 ± 8(stat) ± M(syst) - i[329 ± b(stat) ± 22{syst)\ MeV. (11) 

Values for both k and _Ko(1430) supersede earlier determinations in Ref. |24j because of the 
improvement in the formulae for the k. For the latter, further systematic errors are included to 
cover uncertainties in the extrapolation to the pole, estimated by changing form factors and a 
possible small contribution from the Krj channel. The n pole has moved from 750I55 — i342 ± 60 
MeV [21] by more than the quoted systematic error on mass because of the improvement in the 
fitting formulae. A systematic discrepancy with LASS data has been cured. There is no apparent 
need to increase the flexibility in the fitting formulae, but it is difficult to estimate systematic 
effects this might have on the k pole. 

The scattering length a for the k component is given by Eq. (5) except for a small correction 
from the sub-threshold contribution of Krj via the denominator of Eq. (4). For the optimum 
fit, the k contributes 0.1860/m,,- to the scattering length, i^" (1430) contributes 0.0086/m^ and 
iTo(1950) contributes 0.0002/rrv, i.e. a total of (0.1950 ± 0.0060)/m^. Weinberg [12] predicts 
0.172/m,,- and corrections from Chiral Perturbation Theory to order p 4 increase this to (0.19 ± 
0.02) /m n [28] . Experiment and prediction are now in good agreement and experiment is now 
better than theory. 

The k pole position also now agrees well with the prediction of Descotes-Genon and Mous- 
sallam [20] from the Roy equations, 658 ± 13 — i(278 ± 12) MeV. They considered only the mass 
range up to 1 GeV and omitted the Krj' amplitude and its sub-threshold contribution. 

For completeness, Table 3 shows additional parameters used in fitting the production data 
of FOCUS and E791. The experimental groups use i^*(980) as a reference amplitude. Table 
3 therefore refers to parameters fitting their tabulated amplitudes for the Kn S-wave. Table 3 
also includes the normalisation factor A norm which scales the BES 2 data for i^o(1430). 



Parameter Value 



A K 


2.98 ±0.09 


A1430 


-0.72 ±0.04 


A1950 


-3.16 ±0.07 




130.9 ±2.3 


01430 


53.9 ±2.6 


01950 


26.6 ±4.1 


Anorm 


0.925 ±0.028 



Table 3: Additional parameters fitted to production data. 



3.3 A check using the dispersive correction m(s) 

The s-dependence of m(s) derived purely from the Kir component of T to tai is shown in Fig. 4, 
after a subtraction at the Kir threshold. The fit to data including m(s) is almost indistinguish- 
able from that of Figs 1 and 2. To check the scattering length, it is necessary to fit m(s) to 
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Figure 4: m(s) evaluated from Eq. (8) with a subtraction at the Kir threshold 



terms in k and k 3 . The derived scattering length is (0.201 ±0.007)/m„-, agreeing with the result 
quoted above within the error. One should note that the dispersion integral for m(s) involves an 
integral from threshold to infinity and is therefore sensitive to assumptions about the behaviour 
of the k amplitude above the available mass range. There could be contributions from further 
inelastic channels such as kg. The amplitude itself 'knows' about such contributions, so the 
determination of the scattering length direct from fitted amplitudes at low mass is likely to be 
the more reliable. 

3.4 How the k amplitude varies off the real s-axis 

The Ktt and tttt amplitudes are analytic functions of s. From the Cauchy-Riemann relations, 



there is a rapid variation of both real and imaginary parts of / with Ims. This leads to a 
rotation of the phase of the amplitude with Im s. On the real s-axis, unitarity requires that the 
phase is zero at threshold, but it moves negative with increasing negative values of Ims. 

This is illustrated in Fig. 5 for four values of Ims. Note the the phase below the pole 
moves from a small value for the full curve at Im m = —0.2 GeV to a large negative value for 
Im m = —0.32 GeV, very close to the pole. This explains the curious result that there can be 
a pole with a real part close to zero at threshold. Oiler [29J has drawn attention to this point 
in a somewhat different way. On the real s-axis, what one sees is similar to the upper part of 
the pole, but rotated in phase. This explains why the phase for real s does not reach 90°. It 
reinforces the fact that M and T of Eqs. (1) and (2) are remote from those of the pole. 



d(Ref)/d(Res) 
d(Imf)/d(Res) 



d(Im f)/d(Im s) 
-d(Ref)/d(Ims) 



(12) 
(13) 
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Figure 5: The phase variation near the pole v. mass, for Im m = -0.2 GeV (full curve), -0.25 
GeV (dashed), -0.30 GeV (dotted) and -0.32 (chain curve). 

4 Discussion of results 

The key features emerging from current data are (i) the k peak near threshold, Fig. 3(b), and 
(ii) the precise form of the phase shift near threshold, Fig. 3(c). The FOCUS collaboration [1] 
illustrate in their Fig. 8 the fact that their measured phases near threshold lie distinctly below 
those obtained from a conventional effective range form. This is of course due to the nearby 
Adler zero, which is absent from the usual effective range form. The dashed curve in Fig. 3(c) 
illustrates the LASS effective range formula, which is more curved near threshold. 

Insight into the nature of o and k is provided by the model of confinement constructed by 
van Beveren and Rupp [30]. In this simple model, confinement is approximated by a harmonic 
oscillator potential (which can be solved algebraically), matched to plane wave states at a 
boundary. The way the model is constructed, it approximately reproduces the effect of the 
Adler zero. With the boundary at 0.7 fm, the model reproduces quite well the parameters of all 
of (980) and /o(980) with a single universal coupling constant. A reasonable phase angle 

is needed between uu and ss to reproduce / (980) and a amplitudes near 1 GeV. The / (980) 
and a (980) are locked to the KK threshold by the sharp cusp in m(s) due to the opening of 
this threshold. The ao(980) does not appear at the r)7r threshold because of the Adler zero close 
to this threshold. It is noteworthy that this model was the first (in 1986) to reproduce the lowest 
scalar nonet [31], with the title: 'A low-lying scalar meson nonet in a unitarised meson model'. 

In this model, a, k, a (980) and /o(980) are a nonet of continuum states coupled to the 
confining potential at its boundary. They are meson-meson states at large r, coupled to qq 
components within the confining potential. Doubtless the model could be improved by adding 
meson exchanges at large r. In more detail, it is also possible that diquark interactions play a 
role via coloured configurations, as proposed by Jaffe [32] . 
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There is a clear analogy between a and k and the weak interaction because the amplitude 
rises linearly with s near threshold. The scale of the electroweak interaction is set by the masses 
of W and Z. If the Higgs boson appears as a broad pole like a and k, dispersive effects due to 
opening of WW, WZ, ZZ and tb thresholds will play an important role [33] . 

I am grateful for discussions with Profs. G. Rupp, E. van Beveren, P. Bicudo and J. Ribiero 
over a period of many years. 
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